Discrete Phase Integral Method for Five- Term Recursion 

Relations 



Anupam Garg* 

Department of Physics and Astronomy, Northwestern University, Evanston, Illinois 60208 

(February 4, 2008) 



Abstract 

A formalism is developed to study certain five-term recursion relations by 
discrete phase integral (or Wentzel-Kramers-Brillouin) methods. Such recur- 
sion relations arise naturally in the study of the Schrodinger equation for 
certain spin Hamiltonians. The conditions for the validity of the phase inte- 
gral approximation are derived. It is shown that in contrast to the three-term 
problem, it is now possible to get a turning points "under the barrier", i.e., in 
the classically forbidden region, as well as inside the classically allowed region. 
Further, no qualitatively new types of turning points arise in recursion rela- 
tions with still higher numbers of terms. The phase integral approximation 
breaks down at the new turning points, requiring new connection formulas, 
which are derived. 

75.10Dg, 03.65.Sq, 36.90+f, 75.45.+J 
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I. INTRODUCTION 



The purpose of this paper is to develop the formalism of the discrete phase integral (DPI), 
or Wentzel-Kramers-Brillouin method, for cases where the recursion relation involves five 
terms. Previous use of this method has, as far as we are aware, been limited to three-term 
recursion relations Jl|-§[. Surprisingly, the extension to five terms is not routine, and entails 
novel physical and mathematical considerations. Further, once this extension is understood, 
no additional concepts are required in dealing with recursion relations involving still more 
terms. 

The physical problem which led the author to consider this extension concerns the mag- 
netic molecular cluster [(tacn) 6 Fe 8 02(OH)i 2 ] 8+ (or just Fes for short). This molecule 
has a total spin J = 10 in its ground state, and crystallizes into a solid where the cluster 
has an approximate D 2 symmetry. The interaction between molecules is very weak, and the 
low temperature spin dynamics of a single molecule are well described by the Hamiltonian 

H = -k 2 J 2 z + (h - k 2 )J 2 x - g/2 B 3 ■ H. (1.1) 

Here, J is a dimensionless spin operator, H is an externally applied magnetic field, and 
ki > k 2 > 0. (Experiments reveal k x ~ 0.33 K, and k 2 m 0.22 K.) 

The spectrum of the Hamiltonian ( |1 . 1|) shows some extremely interesting features as a 
function of the applied field H. In particular, one finds a large number of diabolical points || 
in the H x -H z plane pj|-|T2|1 , which have also been seen experimentally JOJ. Exactly as in the 
spectrum of a particle confined to a triangular box , some of the diabolical points (those 



arising when H||x or H||z) can be related to a geometrical symmetry but others can 



not. The problem was first studied [10| by instanton methods when H||x, but this method 
is much harder to apply for general field orientations, and the phase integral method proves 
to be simpler. The existence of diabolical points turns out to depend critically on having 
five terms in the recursion relation that we shall describe shortly, and three terms would 
never lead to such points. The calculations which pertain specifically to Fe 8 are described 



elsewhere [15 , but it seems worthwhile to present the formal aspects of the work separately, 



as they are more generally applicable. 



A. Heuristic discussion of the DPI approximation 

It is useful to continue with the above example in order to introduce the DPI method. 
The starting point of the procedure is to write Schrodinger's equation in the J z basis. Let 
= E\ip), J z \m) = m\m), (m\ip) = C m , and (m\H\m') = t m ^ m i . Then, 

m+2 

^ t tm,nC n EC m . (1.2) 
n=m— 2 

The diagonal terms (t m ,m) in the above equation arise from the J z and J Z H Z parts of 7i, 
those off-diagonal by one (t m , m ±i) from the J X H X and J y H y parts, and those off-diagonal by 
two (t mtm ± 2 ) from the J x part. 

The DPI method is applicable to a recursion relation such as ( |1.2| ) whenever the t m ^ m ± a 
(a = 0, 1, 2) vary sufficiently slowly with m. A physical analogy may be made with an 
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electron hopping on a lattice with on-site energies t m>m and nearest-neighbor and next- 
nearest-neighbor hopping terms £ m , m ±i and t m , m ±2- If these quantities were independent of 
m, the solutions to Eq. ( |1.2| ) would be Bloch waves C m = exp(igm), with an energy 

E = w m + 2t m>m+ i cosg + 2t m , m+2 cos2g = E(q), (1.3) 

where we have written w m = t m>m to highlight the physically different role of the on-site 
energy from the hopping terms. We shall use the notations w m and £ m , m interchangably. If for 
fixed a, the t m ^ m+a vary slowly with m (where the meaning of this term remains to be made 
precise), we expect it to be a good approximation to introduce a local Bloch wavevector, 
q(m), and write C m as an exponential e 1 *, whose phase $ accumulates approximately as 
the integral of q(m) with increasing m, in exactly the same way that in the continuum 
quasiclassical method in one dimension, one writes the wavefunction as exp(iS(x)/h), and 
approximates S(x) as the integral of the local momentum p(x). 

It is obvious that the above approximation will entail the replacement of various sums 
by integrals, and to that end, we introduce smooth functions t a (m) of a continuous variable 
m as extensions of t m ^ m+a such that whenever m is an integer [ |T6|j , 

t a {m) = (t m , m+Q + t m , m _ a )/2, a = 0,1, 2. (1.4) 

We will try and choose these functions so that their derivatives are small. The precise 
way in which this is to be done will be discussed later, but supposing that we have been 
successful in finding such functions, we can seek to approximate C m in exact parallel with 
the continuum phase integral approach. The form of the solution that emerges, and which 
readers will readily appreciate from knowledge of the continuum case, is given by 

C m ~ exp (i I q{m')dm'\ , (1.5) 

Jv(m) ^ J ' 

where q(m) and v(m) obey the equations 

E = w(m) + 2t\(m) cosg + 2t 2 (?7i) cos(2g) = 7i sc (g, m), (1.6) 
v{m) = dH sc /dq = —2 sin q(m)(ti(m) + 4t 2 (m) cosg(m)). (1.7) 

[Just as for the matrix elements, we define w(m) = t (m), and use the notation w(m) when 
we want to emphasize it is as an on-site energy] The interpretation of these equations is 
exactly the same as in the continuum case. Thus, Ti, sc (q, m) is a semiclassical Hamiltonian, 
q(m) is a local wavevector as already mentioned, and v{m) is the associated semiclassical 
electron velocity. We shall refer to Eq. ( |1.5|) as the basic DPI form. Equation ( |1.6| ) is 
the eikonal or Hamilton- Jacobi equation, while Eq. ( |1.7|) is the discrete counterpart of the 
transport equation. The presence of the lattice shows up in the q dependence of H sc (q,m) 
through periodic functions, whereas in the continuum case, such dependence is typically of 
the form q 2 . 

As discussed by Braun |fJQ , the DPI approximation has been employed in many problems 
in quantum mechanics where the Schrodinger equation turns into a three-term recursion 
relation in a suitable basis. All the types of problems as in the continuum case in can then 
be treated — Bohr-Sommerfeld quantization, barrier penetration, tunnel in symmetric double 
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wells, etc. In addition, one can also use the method to give asymptotic solutions for various 
recursion relations of mathematical physics, such as those for the Mathieu equation, Hermite 
polynomials, Bessel functions, and so on. The general procedures are well known and simple 
to state. For any E, one solves the Hamilton- Jacobi and transport equations to obtain q(m) 
and v(m), and writes C m as a linear combinations of the independent solutions that result. 
The interesting features all arise from a single fact — that the DPI approximation breaks 
down at the so-called turning points. These are points where v (m) vanishes. One must 
relate the DPI solutions on opposite sides of the turning point by connection formulas, and 
the solution of all the various types of problems mentioned above depends on judicious use 
of these formulas. 

In this paper we will extend these ideas to five-term recursion relations, focussing espe- 
cially on those features which arise over and above the three-term problem. Now, for any 
given E, there will be four DPI solutions ( |1.5|) , while in the three-term case there are only 
two, because the Hamilton- Jacobi equation ( |1.6| ) is a quartic in e iq . These solutions will also 
break down at turning points — points where v(m) vanishes. In contrast to the three-term 
case, we shall see that there are new types of turning points. It is these turning points which 
are responsible for the diabolical points in the spectrum of Fes- The three-term problem 
turning points are analogous to those in the continuum quasiclassical method, but the new 
ones that we will find are not. In fact, they can only be described as lying "under the 
barrier" from the continuum viewpoint. These new or irregular turning points require new 
connection formulas, which it is our goal to provide. 

The plan of our paper is as follows. We will examine the DPI approximation carefully 
in Sec. [TT], and see how it fails when v(m) vanishes. The precise width of the failure zone is 
discussed in an Appendix. We will examine these failure or turning points in Sec. Ill, and 
see how the concept must be extended beyond the three-term and continuum cases. We will 
find that a turning point need not be a limit of the classically allowed motion, and we will 
categorize the different types of turning points that arise. We will conclude in Sec. [TV] by 
deriving connection formulas at the new turning points. 

We will limit ourselves to problems where the matrix t m n is real and symmetric, t m>n = 
t nim , as it simplifies the analysis, and yet suffices to bring out all essential physical points. 
In our spin example, this means that we only consider fields in the x-z plane. The extension 
to complex Hermitean matrices is cumbersome to carry through, but presents no difficulty 
of principle. We shall continue to couch our discussion in quantum mechanical language, 
thinking of E as an energy eigenvalue, although from the mathematical viewpoint, this is 
not strictly necessary. 



II. THE BASIC DPI APPROXIMATION 

In this section, we will examine the DPI approximation in more detail. The argument 
proceeds in close analogy with the continuum case. We begin by restating the approximation 
in a slightly different way Pjjl2"|. Dividing Eq. ( |1.2| ) by C m , and writing Cm+i = Cm+i/Cm, 
we get 

trn,m—2Cm—lCm ^m,m~lCm ^m,m ^rn,m+lCm+l tm,m+2Crn+lCin+2 E. (2.1) 
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If £m,m+a for fixed a is almost the same over some range K ^> 1 of m's, then we will get 
almost the same numerical equation for the £ m 's over this range, and we will clearly obtain 
a good approximate solution if we replace ratios like Cm+i/Cm by unity. This leads to 

which is a solvable quartic equation in ( m (which is exactly like the factor e tq( ^ introduced 
in Sec. |). The corresponding approximation for C m is 

m 

C m « I] C*, (2-3) 

k=m a 

where m a is a suitable starting value. One could now use this approximation to find the 
ratio Cm+i/Cm, substitute this value for the ratio in Eq. fl2.1|) , and solve again for ( m . The 
process could be further iterated if desired. 

The product in Eq. ( |2.3| ) is more easily evaluated as a sum by taking logarithms. Further, 
if the £ m 's do not vary rapidly from one m to the next, the sum will be well approximated 
by an integral. This makes it necessary to introduce continuum extensions of the matrix 
elements t mim+a . We turn therefore to this problem, and discuss the condition for slow 
variation more clearly. We are seeking functions t a (m) such that [16] 

t a (m) = {t mim+a + tm, m -a)/2, (« = 0, 1, 2), (2.4) 

whenever m is an integer. There are infinitely many such functions, and we restrict them 
by imposing further conditions on a certain number of their higher derivatives. For ti(m), 
e.g., we could also demand (using dots to denote derivatives with respect to m), 

t\{jn) tm,m+l ^"m,m—li (2-5) 
ti(l7l) = | (t m +l jm +2 — tm,m+l ~~ tm,m-l + ^m+l,m+2) • (2-6) 

Similar conditions can be imposed on to( m ) an d ti(m). In general, we need conditions only 
up to some small order for practical applications. Up to the degree of approximation in 
Eq. ( |1.5| ), e.g., we need only go up to second derivatives. 

Since the matrix elements are assumed to vary slowly, we want these derivatives to be 
small, and this condition is best codified in terms of a small parameter e, which plays the 
same role as h in the continuum case, such that t a is formally of order e°, t a of order e, t a 
of order e 2 , and so on. For spin Hamiltonians such as Eq. ( |Llp , this parameter is 1/J. The 
quasiclassicality conditions then read 

dta _ q ( t a( m ) \ d2t * _ Q f t »( m ) \ ( 27 ) 
dm \ J ) dm 2 \ J 2 ) ' 

etc. We will continue to use 1/ J as a generic small parameter in the rest of our analysis. A 
problem in which one can not find functions t a (m) obeying Eq. ( |2.7| ) will not be amenable 
to a phase integral approximation. 

One small point should be kept in mind while judging orders of smallness in the spin 
problem, and others like it. There are natural algebraic expressions for the to,(m)'s in which 
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m appears only in the combination mj J or m/[J(J + l)] 1 ^ 2 - Thus, derivatives with respect 
to m are automatically smaller by an order J" 1 . As J —>■ oo, however, the classical quantity 
is not m, but raj J. Thus the ratio mj J should be regarded as a quantity of order unity and 
not J" 1 . 

With this lengthy preamble, we are ready to solve our recursion relation. As in the 
continuum case, we make the exponential substitution, 



C „ 



(2.8) 



The first approximation is obtained if we assume that $(m) varies slowly so that $(m) may 
be neglected. Then C m ± a ~ C m e ia ® , and substituting this into the recurrence relation ( |1.2|) 
with t m ^ m ± a ~ t a (m), we obtain the Hamilton- Jacobi equation ( |1.6|) with $(m) = q{m). 

To proceed more systematically, we look for a solution for $(m) as a series in inverse 
powers of J: 



$ = $0 + $1 + $2 + 



(2.9) 



where, 



0(J l ~ n ), 

0($ n /J), 

0($„/J 2 ). 



(2.10) 



and so on. The successive inverse powers of J in the derivatives are expected since we 
expect the $ n to be simple functions of the t m ^ m+a , and we shall soon see whether or not 
this expectation is fulfilled. $o is our zeroth order approximation above, and so we set 



$n(m) 



q(m) 



(2-11) 



from the outset. As in the continuum case, we need to keep terms up to $2 m order to 
decide if the approximation is succeeding . 



Up to terms of order 1/ J 2 relative to the leading one, we have 
C m ±i = C m exp 



±q ± $1 ± $ 2 + ~(q + $1) ± \q + ■ 
2 



1 + 



i(g±2<i> 1 ) 2 ±z$2 + ^i±^ + 



(2-12) 



We now wish to substitute this form into our recursion relation. We would like to use the 
continuum forms t a {m) instead of the discrete matrix elements t mjm+a . Just as there are 
infinitely many continuous functions which we could take, there are many approximants for 
t m ,m+a- Note in particular, that there is no unique way to "solve" Eqs. ( |2.4] )- (|2.6|) for the 
matrix elements in terms of the continuous functions. The simplest procedure is to take 
tm,m = t (m), t m>rn ±i = ti{m ± |), and t m>m±2 = t 2 (m ± 1). For £ m , m± i, in particular, a 
Taylor expansion of this approximation gives 



t m ,m±i = ti(m) ± ^i(m) + \ 'ti(m) + 
2 o 



(2.13) 
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where the error is of order J 3 . Therefore, 



where 



tm, n C n = ti(m)ai(m) + ii(m)a 2 {m) + ii(m)a 3 (m) + 

n=m±l 



aAm) = 2C„ 



cos q + - (q cos q + 2i<&\ sin g) (q 2 + 4$ 2 — 4i$i) cos q 



— ( 3g$i — 6i$2 — iq ) sin g 



a 2 ( m ) = iC) 
, , 1 



sin q + -g sin g + $x cos g 



C m cosg. 



(2.14) 



(2.15) 



Note that we have kept only terms up to order 1/J 2 in ckj, 1/J in a 2 ; an d m a 3i since 
relative to t 1 (m), i 1 and £i are of order 1/J and 1/J 2 , respectively. 
Similarly, 



WCn = t 2 {m)^ 1 {m) + i 2 (m)f3 2 {m) + i 2 (m)f3 3 {m) + 

n=m±2 



(2.16) 



where 



A(m) = 2C n 



cos 2g + 2z (g cos 2g + z4> x sin 2g) - 2 (g 2 + $ 2 - cos 2g 



2i / ■ • x 
— — (Gq§i - 3i$ 2 - 2igj sin2g 

I3 2 {m) = 2iC m sin 2g + 2zg sin 2g + 2$i cos 2g 
(3 3 (m) = C m cos2g. 



(2.17) 



We now substitute these relations into the recurrence relation (|1.2|), and equate equal 
powers of J. The terms of order J° obviously give Eq. ( |1.6| ), while those of order J -1 give, 
after some work, 



— 2$! (t\ sin q + 2t 2 sin 2g) = — %— — (t\ sin g + 2t 2 sin 2g) 



(2.18) 



The left hand side of this equation is &iv(m) [see Eq. ( |1.7|) 1, while the right hand side may be 
written as (i/2)dv(m)/dm. Thus, integration yields (after a suitable choice of an indefinite 
constant) 



$i(m) 



lnv(m) 



(2.19) 



Note that as assumed in Eq. ( 2.1C ), $i = 0(J°). Pending a demonstration that Q 2 (m) is 
negligible, we have arrived at the basic DPI form ( |1.5[ ), which we now see as the first two 
terms of an asymptotic expansion of $ in inverse powers of J. 

The equation for Q 2 (m) is considerably more involved. After some analysis, we find 
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where we have defined 



d$2 1 d 2 r i d$i 

£ = |- - r t (2.20) 

c?m 8 dm 2 4 c?m ' 



* * 1 /" ti + 16*2 cos g .. , /ooin 

$ 2 = $2 H / -qdm, (2.21) 

24 i *! + 4* 2 cosg ' V ; 

/ / \ \ h cos q + 4t 2 cos 2q 

rim.qim)) — . 2.22) 

v Hy " ti sin q + 2*2 sin 2g v ; 



Thus, 



$ 2 = f ^— — ^ 2 CQS ^ qdm H / r-^—- \nvdm. (2.23) 

24 i *i+4* 2 cosg 8 dm 8 J dm 2 K J 

Based on power counting, this is indeed of order J -1 , since q and d 2 In vim) j dm 2 are of 
order J~ 2 . Thus, |$2| $i = 0(J°). However, it is plain that this condition is violated 
whenever v(m) approaches zero, for then both lnv(m) and r(m) diverge. To find the actual 
magnitude of $2> we need to know how q(m) and v(m) behave near a turning point. This 
behavior is found in the next section [see Eqs. (|3.14j ) and ( |3.15| )1. The magnitude of <3> 2 is 



estimated in Appendix [A], where we show that the width of the zone where DPI fails is of 
order J 1 / 3 . 



III. TURNING POINTS 

We now turn to a study of the points where v{m) = &H sc (q,m)/dq = 0. We shall call 
all such points turning points in analogy with the continuum case. In contrast to that case, 
however, we will find that turning points are not just the limits of the classical motion for 
a given energy, once the notion of the classically accessible region is suitably understood. 

Since we must also obey the eikonal equation ( |L6p in addition to the condition v(m) = 0, 
at a turning point both m and q are determined if E is given. Setting v = in Eq. (|L7D, 
we see that we must have either q = 0, ot q = tc, or q = q*(m), where 

cosg*(m) = — *i(m)/4* 2 (m). (3.1) 

Substituting these value of q in the eikonal equation, we see that a turning point arises 
whenever 



E = U (m), U w (m), or C/*(m), (3.2) 



where, 



U (m) = H sc (0, m) = w(m) + 2ti(m) + 2* 2 (m), (3.3) 

U v (m) =n sc (n,m) = w(m) - 2*i(m) + 2* 2 (m), (3.4) 

t 2 (m) 

U*(m) = H sc (q*,m) = w(m) - 2* 2 (m) - \ (3.5) 

4* 2 (m) 

Note that q*(m) may be complex for some m, but since cosg* is always real, U* is real for 
all m. We shall refer to these three energy curves as critical curves. 



8 



In the one-dimensional continuum case where 7i sc = (p 2 /2m) + V(x), the condition 
v(x) = is equivalent to E — V(x). The latter condition marks the edge of the classically 
allowed region, E < V(x). Let us recall why this is said to be so. A particle at a point x 
where V(x) > E must be ascribed an imaginary momentum. To understand the analogous 
condition in the discrete case, let us return to the analogy of an electron in a one-dimensional 
lattice with m-independent matrix elements t m , m+a . Equation ( |1.3|) gives the dispersion 
relation for an energy band E(q). The classicaly allowed range of energies is now defined 
by the limits of this band since (provided one is not too close to a band edge) a spatially 
localized electron with a mean energy in this range can be constructed as a wavepacket out 
of Bloch states with only real wavevectors. Alternatively, we could say that solutions with 
energy in the allowed range correspond to travelling waves, while solutions outside this range 
correspond to evanescent waves. These notions continue to be valid when t a (m) are slowly 
varying with m. The dispersion relation may be taken as 7i sc (q, m) for any fixed value of m. 
The band is now m-dependent, and so in particular are the band edges, which we denote by 
U-(m) (lower edge) and U + {m) (upper edge). 

To find the band edges, we note that by definition, the energy at an edge is either a 
minimum or maximum, so &H sc (q, m)/dq = 0, i.e., v(m) = 0. [The converse is not true, i.e., 
v(m) = need not always define a band edge.] Thus the wavevectors at the band edges are 
again to be found in the set q = 0, 7T, and q*, with the answers depending on the signs and 
magnitudes of t\ and t 2 . [We do not distinguish between q = it and q = —tt, or between 
q = q* and q = — q*, as H. sc (q,m) is an even function of q] To narrow down the number of 
cases to be considered, we observe that the gauge transformation C rn — > (— l) m C m changes 
the sign of t±. Hence, we may assume t\ < without loss of generality. Consider the case 
t 2 > first. Then the upper band edge is always located at q = tt, while the lower band 
edge is located at q = if ti/4t 2 < — 1, and at = cos _1 (— ti/4t 2 ) if ~ 1 < ti/^t 2 < 0. The 
case t 2 < is completely analogous. Now the lower band edge is always at q — 0, while the 
upper band edge is at q = tt for £i/4i 2 > 1, and at q* = cos _1 (— £i/4i 2 ) for < ti/4t 2 < 1. 
The different types of band energy curves that can arise are illustrated in Fig. |l|. In the rest 
of this paper we shall carry out the analysis in detail assuming t 2 > 0. The case t 2 < is 
easily treated in parallel, and we shall only give the final results where these are significantly 
different. 

The conditions Eq. ( |3.2| ) for a turning point are analogous to the requirement that 
E = V(x) in the continuum case. Correspondingly, it helps in visualization to draw all 
three critical curves versus m, and a horizontal line indicating the energy. Any intersection 
of this line with a critical curve is a turning point. See, e.g., Figs. and |3|. 

Not every critical curve need be a band edge curve, however. Since we have chosen 
t\{m) < and t 2 {m) > 0, the upper edge curve U+(m) is the same as U n (m) for all m, 
but the lower edge curve, U-(m), may be Uo(m) for some values of m, and U*(m) for other 
values as discussed above. It turns out to be useful to introduce a dual labelling scheme for 
the critical curves and write 



We have already noted that U n (m) = U + {m). The subscripts i and / stand for "internal" and 
"forbidden", since in the first case above, Uo(m) lies inside the classically allowed energy 



Uo(m) = Ui(m), U*(m) = t/_(m), if g* G (0, tt), 
U (m) = U-(m), U m (m) = U f (m), if g* ^ (0,7r). 



(3.6) 
(3.7) 
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range, while in the second case, U*{m) lies outside this range. The turning points with 
E = Ui and E = Uf have no analogues in continuum quantum mechanical problems. 

Before turning to a classification of the various turning points, however, it is useful to 
record some further properties of the critical curves. The first property is that Uo(m) > 
Z7*(m), since 

U (m) - U*(m) = —L-it^m) + 4t 2 (m)) 2 . (3.8) 
4*2 {m) 

Differentiating this equation with respect to m, it follows that the case of equality, U (m) = 
U*(m), happens at a point where both curves have a common tangent. These facts are 
illustrated in Fig. f|. Further, at the point of contact, which we denote by m*, t\(m) / Atzirn) = 
— 1, which is precisely the condition derived above for the lower band edge to change from 
q = to q = q*. 

The second property provides an alternative way of viewing the condition E = U*(m). 
Solving the eikonal equation Q1.6Q for cos q we obtain 



cos,H = ^'"'^y^ 1 " 1 ; (3.9) 

f(m) =w{m) -2t 2 {m) - E. (3.10) 

Since cosg = — £i/4£ 2 at q = q*, the discriminant in Eq. fl3~9[) must vanish, and we must 
have 

t\{m) =4t 2 (m)(w(m) - 2t 2 (m) - E) (q = q m ). (3.11) 

It is easily verified that this equality is identical to E — U^m). 
We now turn to discussing the different types of turning points: 

Type A: E = U-{m) when LL = Uq. See, e.g., Fig. ||. The region m < m c is classically 
allowed. This is analogous to what happens in the conventional continuum quasiclassical 
method — the turning point is located at the boundary of the classically accessible region 
for the energy given. For m just less than m c , there are two solutions of the Hamilton- 
Jacobi equation ( |1.6| ) with q w [(E — U-(m)) / 'a\ 1 / 2 , where a = —(ti + 4t 2 ). For m just 



greater than m c , these values of q continue on to the imaginary axis. The corresponding 
wavefunctions C m change from slowly oscillatory for m < m c to exponentially growing and 
decaying for m > m c . The connection formulas for these solutions are exactly like those in 
the continuum case, and may be derived as in Refs. Note that the other two solutions 

of the Hamilton- Jacobi equation evolve smoothly, and the corresponding DPI wavefunctions 
C m do not need to be "connected" across this turning point. 

Type A: E = U + (m). See Fig. ||] again. The region m > rrib is classically allowed. This 
case is physically very similar to type A in that the turning point is at the boundary of 
the classically allowed and forbidden regions. Now, however, q w 7r in the transition zone, 
so the wavefunctions contain a rapidly oscillating factor (— l) m in addition to all the other 
variation. Although the connection formulas can be derived from those for type A turning 
points by means of the transformation C m — > (— l) m C m , as shown in Ref. 0, their detailed 
form has a very different superficial look. 



10 



Type A 1 : E = Ui(m). Consider Fig. 0, and the energy E shown there, which intersects 
Ui(m) at m — m c . In Fig. |5| we sketch energy bands for this problem for several values of m. 
For m just less than m c , H sc (q, m) = E in just two places, which we denote by ±q(m). (We 
do not show the solution —q(m) explicitly.) For m just greater than m c , two new intersections 
develop at ±g'(m) ~ 0. Thus for m > m c , our wavefunction consists of a sum of four basic 
solutions ( |1.5|) , all oscillatory, while for m < m c we have two oscillatory solutions [associated 
with ±q(m)], and two exponentialy decaying or growing solutions [associated with ±q'(m)]. 
The latter solutions must be related across the turning point by connection formulas, which 
are completely identical to those for type A. The fact that we may have exponentially 
decaying and growing solutions inside a classically allowed region is unexpected from prior 
experience with the continuum quantum mechanical problems, and underscores the point 
that this turning point has no analogue there. We are unaware if it has ever been considered 
in other physical situations where a continuum phase integral approach may be applied. 

Type B: E = Uf(m). This turning point is perhaps the most interesting of all. Since 
the energy now lies outside the classically allowed range, the point lies "under the barrier" , 
and at it, q must be purely imaginary. Consider Fig. [|, and the energy E, which intersects 
Uf at m — m c . We see that for m < m c , there are two solutions to the Hamilton- Jacobi 
equation ( |1.6| ) with q = ±m, where k is real. For m > m c , these solutions acquire a real 
part as well, so that C m [See Eq. ( |1.5|) 1 changes from a decaying (or growing) exponential 
to an oscillating solution with an exponentially decaying (or growing) envelope. As for type 
A' points, this behavior has no analogue in continuum quantum mechanical problems, and 
we are unaware of prior analyses in other contexts. We shall derive connection formulas for 
this case in Sec. |TV| . 

Type B'\ E = C/_(m) when C7_(m) = U*(m). This turning point is like type A in that the 
energy lies at the lower limit of the classically allowed region, but like type B in that q ^ 0. 
The solutions which must be connected are purely oscillatory on the classically allowed side, 
with q ~ ±5*, and oscillatory exponentials (growing or decaying) on the forbidden side. The 
connection formulas are similar to those for case B. 

Our nomenclature for the turning points may have become evident to the reader. The 
letter A indicates that q = or n, while B indicates that q = q* . A bar designates cases 
where q has a value close to 7T, leading to a oscillatory factor in C m close to (— l) m , and a 
prime indicates cases where U* is either U~ or U+ at the turning point value of m. Thus in 
the case t% < 0, t 2 < 0, we would have turning points of type A, A, A', B, and B' . Note 
that the mathematical aspects of the turning point, i.e., connection formulas, are governed 
by the value of q at the turning point, but its physical nature is governed by whether the 
energy lies at the boundary (cases A, A, B', and B'), in the interior (cases A' and A'), or in 
the exterior (cases B and B) of the classically allowed range. 

It will now be apparent that in problems with three-term recursion relations, where 
7i sc (g,m) = w(m) + 2ti(m) cosg, and v{m) = —2ti(m) sing(m), the only turning points are 
at band edges, with q = or n, i.e., of type A or A'. It is also apparent that no new points 
are involved in problems with further neighbor hopping, i.e., recursion relations with seven 
or more terms. Turning points are encountered whenever the energy lies on a critical curve, 
where this term now describes all curves in the E-m plane on which &H sc (q,m)/dq = 0. 
As m is varied through each turning point, two roots of the Hamilton- Jacobi equation for q 
approach each other parallel to either the real or imaginary axis, coalesce, and move apart 
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in the orthogonal direction |17 |. Consider for example the band structure in Fig. |6|, which 
could arise from a recursion relation with seven or more terms. The accompanying critical 
curves are shown in Fig. [7]. For the energy E shown, if m < m a , there is only one real 
solution for q, i.e. q + . (We do not explicitly mention negative values of q.) At m = m a , a 
new real value of q, q2, enters the picture, and splits into two solutions qi- and q2+ as m 
increases further. 

We conclude this section by finding the general behavior of q(m) and v(m) near a turning 
point to = m c . In its vicinity we may write 



H sc (q,m) = H sc (q,m c ) 



on. 



to 



dm 



+ 



(3.12) 



= 0, and so H sc (q,m) 
expand the the right hand side of Eq 



By definition, however, v{m c ) - 
vanish. If we write q{m c ) = q c , 

q — q c as well, and retain only the leading non vanishing terms in q — q c and to — m, 
obtain 



E, and d[H sc (q,m c ) — E]/dq both 
T2|) in powers of 
wc 



H sc {q, to) - E ps a(q- q c ) 2 + b(m - m c ) + 



(3.13) 



where a and h are constants. If we regard w(m) and t a (m) as being of order J°, then by 
Eq. ( |2.7|) , we have a = 0(J°), and b = 0(1/ J), and therefore 



q(m) - q c ~ [(m - m c )/J] 1/2 , 
v (m) ~ [(to — to c ) / J] . 



(3.14) 
(3.15) 



These formulas prove useful when connection formulas are derived. Their usefulness is 
limited, however, if there is another turning point very close to m c . In this case we should 
keep terms of order (to — to c ) 2 . The requisite analysis is very similar to that of quadratic 
turning points [IS] in the continuum case, but we shall not have any occasion to pursue it 
further. In all the calculations we have done for the Fes or other spin problems [|B5 , we have 
been able to sidestep the associated quadratic connection formulas by directly matching the 
solutions in the forbidden region to solutions of the Schrodinger equation for a harmonic 
oscillator. 



IV. CONNECTION FORMULAS FOR FORBIDDEN REGION TURNING 

POINTS 

We turn at last to the problem of finding connection formulas at the turning points. The 
formulas for points of type A, A, etc. are quoted by Braun 0, so we will only consider points 
of type B and B. Our procedure is a small modification of that used by Schulten and Gordon 
0. Suppose the turning point is at to = to c , and q(m c ) = q c . As shown in Appendix |A], the 
DPI solution ( |1.5| ) fails in a window Ato = \m — m c \ < 0( J 1 / 3 ), which we shall refer to as the 
failure zone. The first step is therefore to find another approximation that holds in the larger 
window (which we refer to as the central zone) Ato J v , where rj > 1/3. To do this we write 
C m as an e tqcm times a slowly varying factor, y m , for which we then derive an approximate 
second order differential equation. The second step is to asymptotically match solutions 
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of this differential equation to the DPI solutions in the overlap zones J n ^> Am ^> J 1 / 3 
on either side of the turning point where both types of solutions are valid. The last step 
is to directly write down the transformation matrix between the coefficients of the linear 
combination of DPI solutions for m < m c to those for m > m c , without having to consider 
the solution in the intermediate zone. We will carry out these three steps only to an order 
necessary to match the solutions to the accuracy represented by Eq. ( |1.5| ), i.e., to order J° 
in the phase $(m) introduced after Eq. (|1.2|) or in Eq. ( |2.8| ). 

Let us assume as before that t x < 0, t 2 > 0, and first consider turning points of type B. 
We will denote quantities evaluated at m = m c by a subscript c: t 1 (m c ) = t lc , ti(m c ) = i lc , 
etc. Let q be pure imaginary for m < m c , and precisely at m c let us write 

q(m c ) = ia 2 Kc, (4.1) 
where k c > and a 2 = ±1. Putting this in Eqs. ( |3.1|) and ( |1.6| ) we have 

ti c = — 4^2c cosh k c . (4.2) 
E = w c + 2t lc cosh k c + 2t 2c cosh 2k c . (4.3) 

To carry out step 1, we write 

C m sa e- a2Kc{m - mc) y{m), (4.4) 

where y(m) -C k c . Assuming that this is so, we write 

C m±k = e-^^-^) e -^ [y ±k y + + . . (4 5) 

This approximation will hold provided y y, and y y throughout the central zone 
\m — m c \ J v . [We anticipate, in fact, that throughout this zone, y(m) ~ J~ 1 y(m) where 
7 > 0.] That this is so and that higher order derivatives can be neglected in Eq. (|4.5| ) will 
be justified post facto. Substituting Eq. ( |4.5|) in Eq. ( |1.2| ) we obtain 



A (m)y(m) + Ai(m)y{m) + A2{m)y{m) w 0, (4.6) 



where 



A 






fe=±l,±2 


A 1 


= E ke~ 




fc=±l,±2 


A 2 


= E 




fc=±l,±2 



6 tm,m+k- (4'7) 

b m,m+ki V 41,0 / 

tm,m+k -E ', (4.9) 

The differential equation that we are seeking for y(m) need only hold in window of width 

0(J V ) with 77 > 1/3 around m c . If we choose r) < 1, we can use the fact that w(m) and 
t a (m) are slowly varying, and use the expansions 

«t fc (m)±-fc* fc (m), (4.10) 

tjfc(m) w t kc + (m - m c )i kc , (4.11) 
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etc. Doing this, we obtain 



where 



A ^ ai , (4.12) 
Ax « a 2 - a 2 (m - m c )b 2 , ( 4 -13) 
A 2 « -<x 2 a 3 + (m - m c )b 3 , (4.14) 



ai = ti c cosh k c + 4t 2c cosh 2k c , (4. 15) 

a 2 = tic cosh k c + 4t 2c cosh 2k c , (4. 16) 

6 2 = 2 sinh K c (tic + 4t 2c cosh k c ), (4-17) 

a 3 = sinh K c (i lc + 4t 2c cosh k c ) = ^b 2 , (4. 18) 

b-s = w c + 2tic cosh k c + 2t 2c cosh 2k c , (4.19) 



and where we used Eq. Q4.3| ) to simplify the expression for A 2 . We can also simplify the 
results for a\ and 63. Using Eq. ( |4.2| ) in Eq. ( |4.15| ), we obtain 



a x = 4t 2c sinh 2 K c > 0. (4.20) 



To simplify 63, we recall that the discriminant in Eq. ( |3.9|) vanishes at m = m c . Expanding 
about m c , we have 

16 

t\(m) - At 2 (m)f(m) = -—a 2 (m - m c )t 2 2c + 0((m - m c )/J) 2 , (4.21) 

where a is a positive constant. Differentiating with respect to m and setting m = m c , we 
obtain 

8 

tictic - 2t 2c f c - 2t 2c / c = - —c?t\ c . (4.22) 

Since / = w — 2t 2 — E, using Eqs. ( |4.3| ) and ( |4.2| ), we obtain f c = 4t 2c cosh 2 K c . Also, 
f c = w c — 2t 2c , and the left hand side of Eq. Q4.22D becomes 

— 2t 2c (w c + 2tic cosh n c + 2t 2c cosh 2k c ), (4.23) 

which by Eq. ( |4.19| ) we recognize as —2t 2c b 3 . Therefore, 

h = 4aH 2c /J. (4.24) 

One can also similarly show that 



a 3 = t 2c sinh k c f ^ 
dm \t 2 



(4.25) 



but this result is not particularly useful. 

Let us now examine the order of magnitude of the various a and h coefficients just 
introduced. We first note that since t ac = 0(J°) and i ac = 0(J _1 ), Eq. ( |4.22| ) implies that 
a = 0(1). It follows that a\ = 0(J°), while a 2 , a 3 , b 2 , and 63 are all of order J -1 . 
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To solve the differential equation ( f4.6|) with the approximations ( |4.12| )-( |4.14j ) for the A's, 
we eliminate the first derivative via the substitution 



y(m) = z(m) exp 



a 2 (m - m c ) - ^cr 2 (m - m c ) 2 



(4.26) 



2ai 

Then z(m) obeys 

aiz(m) + 63(777. — m' c )z m = 0, (4.27) 
where we have dropped a term of order (m — m c ) 2 / J 2 in the coefficient of z m , and where 

63 = 63 + ^, (4.28) 
a 2 

m > =m c + — 2 —. (4.29) 
4ai6 3 



This is Airy's differential equation, and the general solution can be written as a linear 
combination of Ai(£') and Bi(£'), where 

(V \ 1/3 

We can now assess the validity of our approximation for y(m). First, let us ask for the 
order of y relative to y. From the known behavior of the Airy functions, this is J _1//3 in the 
failure zone Am < J 1 / 3 , and of order (Am/ J) 1 / 2 in the overlap zone J v ^> Am ^> J 1 / 3 . 
We thus see that the higher order terms in Eqs. (|4.12|) - (|4.14|) are indeed smaller than those 
retained. In the same way the terms dropped in the differential equation for z(m) can be 
seen to be small. Second, the d 3 y/dm 3 term which was neglected in the differential equation 
( |4.6| ) is of order (Am) 3//2 / J 1 / 2 in the overlap zone. Since 77 < 1, this term is smaller than 
the y term, and higher order derivatives are smaller still. 

The next step is to match the solution (|4.4|) and ( |4.26j ) using the known asymptotic 
forms of the Airy functions Ai and Bi, onto the DPI forms for m < m c and m > m c . 
The matching zones can be taken to be any regions in which J v ^> Am ^> J 1 / 3 , where 
1 > T] > 1/3. The leading behavior of the Airy functions is either an exponential or a sine 
or cosine of (Am) 3 / 2 / J 1 / 2 . We can ensure that the term ^(m — m c ) 2 /Aa\ in the exponential 
in Eq. (|4.26|) is inconsequential on both sides if we choose Am J 1 / 2 . Accordingly we take 
the matching zone as 

J 1/2 > \m-m c \ > J 1/3 , (4.31) 

In this zone we can approximate y(m) w z(m), and also neglect the small corrections in 
Eqs. ( |4.28| ) and fl4.29| ). This amounts to saying that 



C m w e~^ m - m ^[ Cl Ai(C) + c 2 Bi(C)], (4.32) 
where c\ and c 2 are two arbitrary constants, and 
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1/3 



C 



(m — m c ). 



(4.33) 



Let us first match Eq. ( |4.32| ) to the DPI solutions for m < m c . In order to treat all four 
solutions at the same time, we rewrite Eq. (|3.9| ) as 



cosgfmj 



-him) + ai[t\{m) - At 2 {m)f{m)} 1 / 2 
At 2 (m) ' 



(m < m c ) 



(4.34) 



where f(m) = w(m) — 2t2{m) — E, and o\ = ±1. The wavevector q(m) thus depends on 
both o"i and a 2 , and hence, so does the velocity v (m). We do not bother to write down these 
forms explicitly, except to note that iaxa 2 v(m) is positive. The DPI solutions may therefore 
be written as 



A 



(71 (72 



2Jiax& 2 v (m) 



expi / q(m!)dm' , 



(4.35) 



where A aia2 is a real constant which is introduced as a notational aid in distinguishing the 
different cases. In the matching zone, it follows from Eqs. (|4.34| , (|4.1|) , ( [4.21| ), and ( |1.7D , 
that 



cosq(m) ~ coshK c + aai((m c — m)/ J) 1 ^ 2 , 

q(m) « ia 2 n c + ioio 2 ^r — ((m c - m)/J) 1/2 , 

smh k c 

v(m) « — 8iao"i<J2 sinh k c |t 2c | (( m c — m)/J) 1 ^ 2 . 



(4.36) 



(We have written \t 2c \ instead of t 2c with a view to including the case t\ < 0, t 2 < 0, which 
we shall consider later.) Therefore, 



.4 



0"lO"2 



4-\/2a)f:2c sinh k c 



m r — m 



J 



-1/4 



cxp 



—<j 2 K c (m — mc) + 



2 a<7i<7 2 (m c — m 



.3/2 



3 sinh rc r 



Jl/2 



Since £ 3> 1 for m < m c in the zone (|4.31 ), we may use the asymptotic forms [21 



Ai(C) 



2^' 



-1/4 



Bi(C) « -=C 



-1/4 



7T 



exp l ^C 3/2 ) • 



(4.37) 

(4.38) 
(4.39) 



Using Eq. Q4.33| ) and comparing with Eq. fl4.37| ), we see that we must use the function Ai 
when <j 2 = — ax, and Bi when a 2 = ax- We can write 



n - KA p -v2K C (m-m c ) / Bi(C), a 2 -ox\ 

l2Ai(C), o- 2 = -ai, 



(4.40) 



with 
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K 



7T 



Ay/2at 2c sinh k c 



r . \ -1/12 

a A j i/4 
03. 



(4.41) 



To connect Eq. ( |4.40| ) with the DPI solutions for m > m c , we make use of the reality 
principle, namely that the solution to the basic recursion relation (|1.2p must be real for 
m > m c if it is real for m < m c . Since the solutions to Eq. fl3.9p now lead to complex q, it 
is clear that we must add two DPI solutions in order to get a real result. To this end we 
define 



9a(?™) = i& 2 K{m) + x{ m )i m > m c 

where «(m) and x( m ) are both real and positive, «(m c ) = /c c , and 

cosh/tcosx = — ti/4|t 2 |, 
sinhKsinx = (4* 2 / - t\) x/2 / A\t 2 \. 

In terms of these quantities, we have 

cosg a (m) = cosh/t(m) cosx(^) — i&2 sinh /■c(m) sin x(m), 

s a (m) = —ia 2 v a [m) = 8^2(^7.) | sinh /t(m) sinx(m) sing a (m) 



(4.42) 



(4.43) 
(4.44) 



(4.45) 
(4.46) 



[As before, we have written |£ 2 (m) | instead of t 2 {m) with a view to subsequently including 
the case £1 < 0, t 2 < 0.] The corresponding quantities for the complex conjugate solution are 
given by qb{m) = q*(m), Sb(m) = s*(m). The DPI solution into which Eq. ( 4.35 ) continues 
for m > m c may therefore be written as 



1 



yjs a (m] 



:exp(i / q a (m!)dm! + iA 

\ J m c 



C.C. 



(4.47) 



where B ai(T2 is real, A fflff2 is a phase to be found, and c.c. denotes "complex conjugate". 

It remains to compare the asymptotic forms of Eqs. ( 4.4C ) and ( 4.47|) in the overlap zone 
( |4.31| ), and thus find A fflff2 and the relation between A aiU2 and B aiU2 . Using the asymptotic 



forms for Ai and Bi, we have 



Cm 



-(T2K c (m—m c ) 



_ C) -l/4 




K-0 3/2 + f 
l(-o 3/2 + f 



0"2 = CTl! 

a 2 = —0i. 



(4.48) 



To write Eq. ([4.47|) in a similar form, we first note that in the matching zone, 



q a (m) w ia 2 n c + 02 



a 



((m-mJ/J) 1 / 2 , 



Sa(m) ~ 8t2 C asinhK 
Therefore, 



sinh k 

f m — m c \ 1 / 2 



V J 



exp 



7T 

2 



cosh/t c ( m — m c \ 1 / 2 



«02— — «02«— . o i T 
sinh K r V J 



i / q a (m')dm' = — cx 2 K, c Am + i— (— C) 3 ^ 2 



(4.49) 
(4.50) 

(4.51) 
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where we have used Eqs. ( f4.20|) , (|4.24j ), and (|4.33| ) to simplify the term in (— C) 3 ^ 2 - Compared 
to this term, the correction proportional to (Am) 1 / 2 in the exponent for s a {m) is of higher 
order in 1/J in the matching zone, and can be ignored. Doing this, and making use of 
Eqs. ( 1Q3D and ( ggg) , we get 



Cm. ~ l=B n 



-0 



— 1/4 — (T2ftcAm 



exp (^("C) 3/2 - ^ + ^A fflff2 ) + c.c. 



(4.52) 



Comparing with Eq. Q4.48|) , we see that B aiU1 = A aiU1 /2, while B, 
A CT1C72 vanishes if U\ - 
summarized as 



A ai -ai- Further, 

•1, and equals ±7r/2 if ±02 = o\ = 1. These results can be 
2-5, 



A 



7T 



(4.53) 
(4.54) 



Putting together Eqs. (|4.35|) , ( |4.47| ), ([4.53|) , and ( |4.54| ), the connection formula may be 
written in the final form 



— = exp 1 



1 2 °"2 ) o : 



q(m')dm' <— C m 

1 / f m n \ 

^=expU / q a (m')dm + + a 1 )a 2 )+ c.c. 



s a (m) 



(4.55) 



The wavevectors q and q a , and the velocity and speed, v{m) and s a , to be used for a given 
set of signs (ci, cr 2 ) are given by Eqs. ( |4.36|) and ( |4.49|) . 

We use the double arrow notation C m — > advocated by Heading |22j and Dingle 
p3| to emphasize the bidirectionality of the connection formula. We refer readers to these 
authors for lucid discussions of this point, but since even as masterly and authoratative a 
text as Landau and Lifshitz |M| states that connection formulas may only be used in one 
direction, it may be worth paraphrasing their remarks briefly. Thus, as stated by Heading, 
the notation means that there is a solution C m to the recursion relation (|1.2[ ) with the stated 
asymptotic behaviors for m < m c and m > m c . And, as stressed by Dingle, the formula 
merely states that a certain exponentially growing on one side of the turning point, which 
is free from the growing component, goes over into a certain oscillatory solution (with an 
exponentially growing envelope) on the other side, and vice versa. Likewise for the growing 
component. It says nothing about whether or not we can use the formula in both directions 
in physical problems where we do not have complete information. If, e.g., a solution with 
o~i = o~2 = — 1 for m < m c contains an admixture of the o\ = —02 = 1 solution of small but 
indeterminable magnitude, then we can say nothing about the amplitude or the phase of 
the oscillatory factor for m > m c . It is in these situations, where information is incomplete, 
that the reservations about the unidirectionality of the connection formulas are relevant. 

The connection formula for the case where t± < and t 2 < can be derived in exact par- 
allel by making minor modifications to the intermediate steps. The final form is sufficiently 
different to be worth giving separately: 
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A 

cr l cr 2 



2\Ji<j\G2V (m) 



expz / q(m')dmf <— C n 

J m r 



1 2 < ^ f7l<T 2 ) ^<J\a.i 



exp (i / q a {m')dm! — — (T^a-A + c.c. 



(4.56) 



>s a (m) 

The quantities v(m) and s a (^) are given by the same formal expressions as before, but now 

q{m) « 7r + ia 2 K c + icri<T 2 ^-; ((m c - m)/J) 1/2 , (4.57) 

smh k c 

g a (m) = 7r + icr 2 K,(m) + x( m )- (4.58) 

In Eq. ( |4.58[ ) , K,(m) and x( m ) are a l so given by the same formal expressions as before, i.e., 
Eqs. (TO ) and (g^§. 
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APPENDIX A: WIDTH OF FAILURE ZONE OF DPI APPROXIMATION 

We have seen in Sec. fl| that the DPI approximation is valid everywhere except near 
turning points. Let us now estimate the size of the region where it fails. Suppose that q = 
or 7r at the turning point, m c . As shown in Eqs. ( |3.14 ) and (|3.15| ), q(m) and v (m) both vary 



as [(m — m c )/ J} 1 / 2 for m near m c . Thus q ~ (m — m c ) 3 ^ 2 J 1//2 , and the first term in $2 is 
of order [(m — m c )JY l/2 . Next, note that r ~ [J/( m — mc)] 1 / 2 , so that the second and third 
terms are both of order J 1//2 (m — m c )~ 3 / 2 . The first term is subdominant to these two, and 
therefore, the condition that $2^1 reduces to 



\m — m c 



> J 1/3 . (Al) 



In the case where the second factor in v(m) vanishes, we have ti + 4t 2 cosg ~ [(m — 
m c )/J] 1//2 , and q — q c ~ [(m — m c )/ J] 1 ^ 2 . Thus the integrand of the first term in $2 varies 
as l/(m — m c ) 2 , and the term itself is of order l/(m — m c ). However, r ~ [J/(m — m c )] 1 ^ 2 as 
before, so that the second and third terms are again of order J x l 2 {m — m c ) -3 / 2 , and much 
greater than the first term. The condition for validity of the DPI form is again given by 
Eq. Q. 

The condition can be written in alternative form by noting that near the region of failure, 
r(m) ~ l/v(m), and that v ~ q. Thus dr/dm ~ q/v 2 , and the last term in $2 can also be 
shown to be of order q/v 2 . Thus the condition for validity can be written as 

q(m) < v 2 {m). (A2) 
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FIGURES 

FIG. 1. Schematic band energy curves arising from five term recursion relations. Curve 1 arises 
if | i i / 4*2 1 > F curve 2 if — 1 < ii/4i<0, and curve 3 if < ti/4i2 < F Fhe dual labelling scheme 
for the critical energies Uo, U n , etc. is shown. 

FIG. 2. Critical energy curves in the case where — ti(m) /4t2(rn) < —1 for all m. Fhe points 
m = nib and m = m c are turning points of type A and A, respectively, for the energy E. 

FIG. 3. Same as figure 2, but with — tx(m) /4t>2(m) 6 [—1,1] for some m. Fhese critical curves 



arise from the Hamiltonian (1.1) for Fes- Note that U* is the lower band edge U- for the central 
m region, and the forbidden critical energy Uf for the outer m regions. Correspondingly, Uq is the 
internal critical energy Ui in the central region, and the lower band edge f/_ in the outer regions. 
Fhe point m = m c is a type A' turning point for the energy E. 

FIG. 4. Magnified view of lower left hand portion of Fig. ^. The points m = mt and m c are 
turning points of type A and B respectively. 



FIG. 5. Band energy curves in the vicinity of a type A' turning point. 
FIG. 6. Possible band energy curve for a recursion relation with seven or more terms. 
FIG. 7. Critical curves for a system with band energy curves as shown in Fig. ^| 
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